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t C"| Abstract 

I Quantum coherence conservation is shown to be achieved by a very high rate of dissipation of an environmental system coupled with a principal 
"£h system. This effect is not in the list of previously-known strategies of noise suppression, such as Zeno effect, dynamical decoupling, quantum 
error correction code, and decoherence free subspace. An analytical solution is found for a simplified model of a single qubit coupled with an 
^ environmental single qubit dissipating rapidly. We also show examples of coherence conservation in a spin-boson linear coupling model with 
^"'a numerical evaluation. 
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' Suppressing decoherence in quantum systems is of growing 
interest in the light of the rapid development of quantum com- 
puting II 1121 . There have been many schemes proposed for this 
purpose in several categories, namely, Zeno effect [3], dynami- 
cal decoupling (or bang-bang control) [4 5 1, quantum error cor- 
rection code (QECC) |6 7 8|, and decoherence free subspace 
(DFS) 0911011 II . Conventional schemes for suppressing deco- 
herence focus on controlling a principal system that is under 
influence of environmental systems. The quantum Zeno effect 
uses a sequence of projective operations (usually projective 
measurements) mapping states onto subspaces of a state space 
|12|. A quantum state is kept inside of a subspace if projec- 
tive operations are applied very frequently before noise kicks 
the state out of the subspace. Dynamical decoupling has been 
studied intensively in the field of nuclear magnetic resonance 
(NMR) from the dawn of the field (see, e.g., Ref. [13]) and 
later extended to other physical systems B4I5I14I15I16II . This 
scheme uses a train of regular short pulses applied to a prin- 
cipal system to cancel time evolutions under noise. The con- 
cept of QECC is to discard a subspace easily affected by noise. 
A QECC uses a space of code words generated from original 
states by adding a certain redundancy so that a recovery from 
corruption is possible. The DFS scheme utilizes a subspace of 
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states that are unchanged by given noise operators. Quantum 
computation is performed in the subspace. 

Studies on coherence conservation so far are mostly based on 
the assumption that we do not have control over environment. 
There may be alternatives available if we focus on the idea of 
controlling environmental parameters. A method to control ef- 
fective coupling using a dynamical control field in the presence 
of time-dependent external field was recently reported by Jirari 
and Potz 1171181 . We pursue the possibility of passive control, 
rather than dynamical control, in a simple picture of decoher- 
ence (See Ref. [19| for conventional decoherence models). 

There are many parameters affecting the decoherence fac- 
tor in realistic models. It is often mentioned that there are the 
strong coupling regime and the weak coupling regime for a sys- 
tem consisting of a principal system and a noise source in gen- 
eral. (A typical case is a qubit under a random telegraphic noise 
|20|.) The two regimes involve significantly different depen- 
dencies of the decoherence factor on the parameters describing 
noise. We need to choose a proper parameter with which one 
can alter the behaviour of the decoherence factor for our pur- 
pose. 

We report in this letter coherence conservation achieved by 
a very high rate of dissipation of an environmental system cou- 
pled with a principal system. A highly dissipative environmen- 
tal system is wiped out before affecting the principal system. A 
clue of the phenomenon was originally found in our numerical 
simulation of bang-bang control of entanglement in a spin-bus 
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Fig. 1. Model for the system consisting of the principal system (system 1) 
and the environmental system (system 2) whose time evolution is governed 
by the Hamiltonian H. System 2 is replaced with a thermal environmental 
system with the dissipation probability p for the time interval r. 

model ||2TI . We investigate the phenomenon in detail by using 
a mathematical analysis and a numerical simulation for simpli- 
fied models. An analytical solution is found for a model of a 
single qubit coupled with an environmental single qubit dissi- 
pating rapidly. The effect is also found in a spin-boson linear 
coupling model by using a numerical evaluation. 

Let us assume that the entire environment is so large that 
the environmental system (system 2) coupled with the principal 
system (system 1) is a part of a large environment and hence 
it is replaced with a thermal environmental system with prob- 
ability p (namely, with some dissipation rate) per certain time 
interval r. Systems 1 and 2 are represented by the density ma- 
trix pi 1 ' 2 ! ; a thermal environmental system is represented by the 
density matrix a. The Hamiltonian affecting the time evolution 
is reduced to the one consisting only of the time-independent 
Hamiltonian H that governs systems 1 and 2 including their 
interaction. This model is illustrated in Fig.Q] For a small time 
interval At, the evolution of the systems 1 and 2 obeys the 
equation 



P ^ 2 \t+At) = e 



-iHAt 



x At p^(t) 



+ (l-x At )Tr 2 p [1 ' 2] (t)®a 



(1) 



iHAt 



where x = (1 — p) x l T and t denotes a certain time step. The 
dissipation rate p can be modified by changing the experimental 
setup under a static control. 

Let us begin with an analytical evaluation of this model in 
the case where the systems 1 and 2 are single-qubit systems. 
For further simplicity, we impose the following conditions. The 
principal system is originally represented by a density matrix 



P [11 (0) 




with < a < 1 and < |6| < y/ a(l — a). The environmental 
system at thermal equilibrium is represented by the maximally- 
mixed density matrix 



P [2] (0) 




The initial state of the total system is set to p^' 2 ^(0) 
/9 [1] (0) <g> p [2] (0). The Hamiltonian H is set to cl z <g> I z 



p^ 2 '(mAi) = 



diag(c/4, -c/4, -c/4,c/4) [here, I z = diag(l/2, -1/2)]. 

With these simplifications, one can find the form of the den- 
sity matrix at time t = mAt (to E {0, 1,2,.. .}) as 

U/2 f m \ 

a/2 g m 

/* (l-o)/2 

\0 g* m (l-a)/2j 

with functions f m and g m depending on to. The functions f m 
and g m obey the system of recurrence formulae: 

/m+l — 77 e ^ [fm "F 9m "F X (fm £/m)] 

{ 

g m +i = -e lcAt/2 [f m + g m - x At (f m - g m )] 

with /o = go = 6/2. This leads to the following recurrence 
formula: 



Krn+2 = (1 + X ) COs(cAt/2)K m+ i - X K m , (2) 

where n m = f m ox g m with / = g = 6/2, } 1 = 6e~ lcAt / 2 /2, 
and g x = 6e lcAt / 2 /2. 

One can derive functions /(t) = liixiAt^o. mAt=tfm and 
g(t) = HniAt^o.mAt=t5m in the following way. By linearliza- 
tion, Eq. (fJJ is put in the form: 

K-m+2 - 2K m+ i + K m - Athlx(K m+ i - K m ) 



+ (At) 2 —n m+1 
+ 0[(At) 3 



(Atf 



(lna;) 2 (K m+ i - n m ) 







Dividing this equation by (At) 2 and taking the limit At 
lead to 







d 2 n(t)/dt 2 - In a; dn(t)/dt + c 2 n(t)/A = 0, 

where n(t) = \iniAt^o,mAt=tK m - The solution of this differ- 
ential equation is 



«(*) = 



+ v K e~ 



with constants u K and v K (k — f or g), and the complex deco- 
herence factor 



lnx± y/(lnx) 2 - c 2 /2 



The real part of r± is nonnegative because In x = ln( 1 — p) /r < 
and | In a; | > |Re-\/ (In a;) 2 — c 2 \. (Here, the square root is 
positive.) 

We need to impose the conditions that k(0) = 6/2 and 
k'(0) — liniAt-»o(^i — Ko)/At. The latter condition can be 
written as — r+u/ — r-Vf = —ibc/4 and —r+u g — r-v g = 
ibc/A. With these conditions, we obtain 

ibc — 26r_ — ibc + 26r+ 



A(r+ -r-Y 
—ibc — 26r_ 



9 ~ 4(r+-r_) 
Consequently, we have 



/(*) = 



ibc — 26r_ 
4(r+ - r_) 



4(r+ -r_) ' 
ibc + 26r+ 
4(r + — r_)' 

— ibc + 2br + 
g 

4(r + - r_) 



<?(*) 



— ibc — 2br. 

~W7 



) 



ibc + 2br + 
4(r+ - r_)' 



One can now write the reduced density matrix of the principal 
system at t as 



^(*) 




with 



-r+t 



Let us investigate r± in details in relation to p. One obvious 
fact is that exponential decay is caused by the real part of 
e~ r±t . Thus we focus on the behaviour of Re r±. We find that 
Re r + increases as — (lnx)/2 when < p < 1 — e~ CT and 
decreases with convergence to zero when 1 — e~ cr < p < 
1. In contrast, Re r_ increases as — (lnx)/2 when < p < 
1 — e~ cr and increases more rapidly when 1 — e~ CT < p < 1. 
This is clearly depicted in Fig. This result suggests that the 
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Fig. 2. Plot of Re r±/c as functions of — (lnx)/c. The decoherence factors 
Re r± increases until they reaches c/2 as the dissipation rate increases. The 
factor Re r+ starts decreasing at the point of — In x = c (i.e., p = l — e~ CT ) 
while the factor Re r_ starts increasing rapidly at this point. 

decoherence factor Rc r + is small for a large dissipation rate 
p > 1 — e~ cr . In addition, Rc e~ r - rapidly converges to for 
p > 1 — e~ cr . Thus the dominant term for such a large p is 
b e~ r+t , which converges to b. This phenomenon can be 
understood physically: the environmental system is wiped out 
quickly before absorbing coherence information of the principal 
system as p approaches to unity. Therefore this effect may be 
called quantum wipe effect. 

In addition, the behaviour of Im r± is rather simple, as shown 
in Fig. [3] They vanish for p > 1 — e~ CT ; this suggests that we 
do not observe oscillation in rj(t) for such a large p. 

As an example, we will see the time evolution of \rj(t)\ for 
r = 1.0 x 10~ 3 s and c = 1.0 x 10 3 Hz. Figure [4] shows clear 
suppression of decoherence for small p and also for large p. 
Suppression of oscillation is found for large p. 

The phenomenon we have seen above can be found also in 
a spin-boson linear coupling model. Let us consider a system 
consisting of a single spin-1/2 (system 1) and a mode of bosons 
(system 2) coupled with the spin system. We assume that the 
mode of the bosons is in resonance with the precession of 



0.5 



0.25 



-0.25 



-0.5 









Im 

Irr 


r + /c - 
r_/c •■ 









































1 2 3 4 5 6 
-(In x)/c 

Fig. 3. Plot of Im r±/c as functions of — (lnx)/c. These imaginary factors 
vanish for — In x > c (i.e., p > 1 — e~ CT ). 
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Fig. 4. Time evolution of \r)(t)\ for several different values of p (0.0, 0.25, 
0.5, 0.75, 0.95, and 1.0) when r = 1.0 X 10" 3 s and c = 1.0 X 10 3 Hz. 

the spin system (the resonance frequency is v). Assuming that 
systems 1 and 2 are surrounded by a bath of bosons of a cavity 
and/or a sample holder, we employ the same model illustrated 
in Fig. Q] 

Let us set the initial state of the spin to 



P [11 (0) = 



The initial state of the bosonic system is (0) = e _/3ff jZ 
where [3 is defined as (3 — (/cbT) -1 (fee is the Boltzmann 
constant and T is temperature), — va)a is the bosonic 
system Hamiltonian (at and a are the creation and annihilation 
operators), and Z is the partition function. The thermal state 
probabilistically replacing is a — p^(0). The initial state 
of the total system is set to p^(0) = p [11 (0) (g) /O [2] (0). Let 
the full Hamiltonian be H = ijW + + H c with the spin- 
system Hamiltonian i/W = vl z and the coupling Hamiltonian 
H c = cl z (a* + a) (here, c is a coupling constant). 

A numerical simulation is used to compute the time evolution 
governed by Eq. (|T). We have taken the parameters: v — 3.4 x 
10 10 Hz, T = l.OmK, c = 1.0 x 10 7 Hz, and r = 1.0 x 10~ 8 s. 
The values of v and T are taken from physically available values 
for an electron spin coupled with an on-resonance bosonic mode 
(a cavity mode and/or a phonon mode of a solid sample) in a 



low-temperature Q-band electron-nuclear double resonance (Q- 
band ENDOR) system. In simulations, At is set to 5.0 x 10~ 10 s. 
The simulation showed a clear suppression of decoherence 

in \(0\pW{t)\l)\ = |(0|Tr 2( o [1 < 2l (i)|l)| for large p as shown in 

Fig.m 




Fig. 5. Numerically computed time evolution of | (0|pW (t)|l)| for several 
different values of p. 

For the next example, we numerically evaluate the decay of 
entanglement of a principal system (system 1) consisting of one 
electron spin and one nuclear spin each of which is coupled 
with an on-resonance bosonic mode (system 2). The measure 
of entanglement that we use here is negativity H22I231 The 
evolution of the systems is governed by Eq. (U. Let us consider 
the following Hamiltonians (the subscripts and 1 denote labels 
assigned to individual spins and corresponding on-resonance 
bosonic modes), (i) The Hamiltonian of the principal system: 
= vqIzq + v\l z \ + AqiI z qI z i with the precession fre- 
quencies vq, v\ and the spin-spin coupling constant A§\ (here, 
I Z i is an I z operator acting on the ith spin), (ii) The Hamilto- 
nian of the bosonic system: = z^oajao + V\a[ax. (iii) The 
Hamiltonian of spin-boson couplings: H c = cqI z q{ci q + ao) + 
cil z \{a\ +a%) with constants cq and c\. The total Hamiltonian 
is H = W + + H c . 

The initial reduced density matrix of the principal system is 
set to 



P [1I (0) 



|00)(00| + |00> (11 1 + |11)(00| + 1 1 1> <1 1 1 



and that of the bosonic system is set to the thermal one p^ (0) = 
a = e-P Hm /Z. We set ^(O) = pW(0) <8 p [2] (0). The 
following constants are used: v Q = 3.4 x 10 10 Hz, V\ — 4.87 x 
10 7 Hz, A i = 1.0 x 10 7 Hz, T = 1.0 X 10- 3 mK, c = ci = 
1.0 x 10 7 Hz, and r = 1.0 x 10~ 8 s. Here, the values of Uq, 
v\, and T are chosen by considering an electron spin and a 
nuclear spin coupled with on-resonance bosonic modes in a 
Q-band ENDOR system. A numerical simulation is performed 
with At = 5.0 x 10~ 10 s. 



1 Negativity Af can be calculated for a density matrix p[ a ' b l of a bipartite 
system consisting of subsystems a and b in the following way: J\f(p ia > b ^) = 
<g) Ax)p[ a ' 6 ] || — l]/2 where Ax is the transposition map acting on the 
subsystem b and hence / ® Ax is a partial transposition map. 



We show the time evolution of negativity of the principal 
system in Fig. [6] An improvement of the entanglement conser- 
vation is found for large p and it is especially conspicuous for 
p > 0.95. 
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Fig. 6. Numerically computed time evolution of negativity of the principal 
system for several values of p. 

In summary, we considered a model in which a principal 
system is coupled with an environmental system that is prob- 
abilistically replaced with a thermal environmental system. In 
such a model, a high dissipation rate of the environmental sys- 
tem results in a conservation of coherence of the principal sys- 
tem. This phenomenon may be called quantum wipe effect. We 
have shown that this phenomenon is found in a simple model of 
qubit-qubit coupling by an analytical calculation. An increase 
of a dissipation rate over a certain threshold value results in co- 
herence conservation while an increase of the dissipation rate 
otherwise results in larger decoherence, as shown in Fig. [2] This 
is reminiscent of the relation between the Zeno effect and the 
Anti-Zeno effect (24), namely that very quick applications of 
projective operations are effective for coherence conservation 
while slow applications of them result in larger decoherence. 
The quantum wipe effect has been also found numerically for 
a model of spin-boson linear coupling. 

The usefulness of the quantum wipe effect is dependent on 
how small the threshold value of the dissipation rate is. A cou- 
pling constant is required to be small to obtain a small threshold 
value. It is thus expected that this effect will be verified exper- 
imentally in future with a weak system-environment coupling. 
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